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We study the spectral statistics for extended yet finite quasi f-d systems which undergo a transi- 
tion from periodicity to disorder. In particular we compute the spectral two-point form factor, and 
the resulting expression depends on the degree of disorder. It interpolates smoothly between the 
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I. INTRODUCTION 

The spectrum of an unbounded periodic system is arranged in continuous bands and the corresponding eigcnfunc- 
tions are extended (unnormalisable). When sufficient disorder is introduced, the system is Anderson localised - the 
spectrum is point like and the eigenfunctions are localised (normalisable) . The transition from a continuous to a point 
spectral measure is a drastic effect, which might have been used to characterize the transition. However, this approach 
is of a limited value, since in practice one always deals with finite systems, where the spectral measure is point-like 
both in the periodic and in the disordered situations. In finite systems, the mean spectral density is independent of 
the degree of disorder. Therefore, for finite systems, the effect of disorder on the energy spectrum can be discerned 
only in the spectral correlations. Indeed, this approach to the characterization of the Anderson transition in three 
dimensional systems was used ||, and the spectral measures were shown to undergo an abrupt change when the 
critical level of disorder is reached. In the present paper we study the spectral statistics for finite quasi 1-d systems 
which undergo a transition from periodicity to disorder. (Quasi 1-d disordered systems of finite length can be either 
"metallic" or "insulating" depending on whether the localization length is larger or smaller than the system length. 
We shall consider only the first case, and the strength of the disorder will be restricted accordingly, to the range 
of values which is sometimes called weak disorder). We shall focus our attention to the spectral two-point form 
factor, and show that it depends very sensitively on the degree of disorder, and derive a universal expression, which 
' interpolates continuously between the periodic and the disordered yet metallic limits. 

The spectral form factor is the main object of our discussion, and it is denned in the following way. The spectrum 
is unfolded by introducing the dimensionless energy e, through the relation de = (d)(E) dE where (d)(E) is the mean 
spectral density. The corresponding dimensionless time r measures time in units of the Heisenberg time in = 2Trh(d). 
We consider a finite spectral interval of length Ae centered at e c , and denote its characteristic function by x(e — e c ). 
Since the mean spectral density of the unfolded spectrum is unity, the number of states in the interval Ae is Af = Ae. 
This energy interval should be sufficiently large so that Af ^> 1 , and sufficiently small so that the mean level density 
and the classical dynamics do not change much as the energy is scanned across it. The oscillatory part of the spectral 
density in this interval is 
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d{e) = x(e- e c ) 
The Fourier transform of this function is 



(1) 



4(r) = / e- 2 " ier d(e) de 

= E^- e c) e " 2 " e,T -^(r) • (2) 
g 

The Fourier transform of the normalized characteristic function is denoted by 6a t (t) and its width is At ~ 1/Ae. 
The form factor is expressed as 

K(r) = ^{\d c (r)\ 2 ) c . (3) 

We use (-) c to denote the spectral average, which is taken over the non overlapping energy intervals located about a 
set of e c values. One can also perform the averaging over any free parameter of the system or over disorder when it is 
introduced. It can be easily shown that (|^) is nothing but the Fourier transform of the spectral two-point correlation 
density S. For a discrete spectrum the normalization in ([}]) is such that the form factor approaches a constant 7 as 
r — > 00 where 7 is the mean spectral degeneracy. 

The expressions for the spectral form factors in the extreme situations of exact periodicity and weak disorder are 
known. In the latter case, when the length of the system does not exceed the localization length, and assuming that 
the Heisenberg time is shorter than the Thouless time the spectral statistics takes the form M 



K(r) = j frV^Tac for r < 1 (4) 
w [ 1 for r > 1 . w 

The factor gx can take the values 1 or 2 depending on whether time reversal invariance is respected or violated, and 
c is the conductivity of the chain. The spectral form factor for periodic systems was recently derived using both field- 
theoretical methods (|] and the semiclassical approximation Since the latter theory is the basis for the approach 
developed in the present paper, we shall describe it briefly to introduce the concepts and the notations which will be 
used in the sequel. 

We consider a chain of N identical chaotic unit cells of length a = 1, with periodic boundary conditions, such that 
the full system shows a discrete translation invariance (Fig. [j](a)). (Alternatively, we could discuss a disordered ring 
configuration which is threaded by an Aharonov-Bohm flux line. This is the system analyzed in . In such a system, 
the classical evolution within a unit cell becomes ergodic after a short time, and one can approximate the classical 
evolution in the entire chain by diffusive evolution. We shall denote the diffusion constant by D. The time it takes 
the diffusive evolution to cover the phase-space uniformly is the Thouless time. 

Due to translation invariance, the quantum spectrum consists of discretised energy bands whose width depends on 
the (dimcnsionless) conductivity per unit cell. It is defined as c\ — 2 / K%{d\) D / 'a 2 , where {d\) is the mean level density 
per unit cell. A few examples of typical bands are shown in Fig. ^[ One can see that for low ci, the bands are flat 
and show little structure. For high values, the bandwidth is of the order of the inter-band spacing, and the bands can 
hardly be recognized if the discretisation is too coarse. 

If the system under discussion is invariant under an anti-unitary symmetry (such as e.g., time-reversal) the bands 
are symmetric about the center and the edges of the Brillouin zone, and the levels are doubly degenerate (7 = 2). 
The reflection symmetry and the degeneracies are broken if the symmetry is lifted, and in this case 7 = 1. 

The quantum spectrum is characterized by two energy scales, the mean intra-band spacing and the mean inter-band 
spacing. The ratio between them is at least N, the number of unit cells. We are interested here in the large N limit, 
and therefore these energy scales are very well separated. Since (d) « (di)N, the spectral correlations which pertain 
to the inter-band scale affect the behavior of the form factor in the range < t < 1/N. The correlations between 
levels in the same band leave their mark on K(t) in the domain 1/N < r < 1. The fact that the spectrum is composed 
of discrete (possibly degenerate) energy levels is expressed in the spectral form factor in the domain 1 < r, where the 
form factor approaches the constant value 7. 

We used different approximations to express the form factor in the three domains mentioned above |Q. 

• 0<t< 1/N : Here one starts from the semiclassical trace formula and employs the "diagonal approximation" 
[0 to write 
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K{t) w g T NTP{r) 



(5) 



The factor is due to the discrete translation symmetry, because of which any generic periodic orbit is replicated 
TV times in the system, g^ stands for the classical degeneracy due to time-reversal (or any other anti-unitary) 
symmetry and it can take the values 1 or 2. P(t) is the classical probability to stay in the same unit cell 
from which the trajectory started, after the time t — tLh ||. Because phase-space is covered diffusively, 

P (t) ~ 1/2 and hence > 



K(t) m g T Ny/N T /2 Cl , (6) 



where c\ is the dimensionless conductivity per unit cell which was introduced above. 

• 1/N < t < 1 : As r increases, the form factor provides information on a finer energy scale. In the vicinity 
of r = 1/N, the energy levels within a single band cannot be resolved, hence K(t w 1/N) takes a value 
which is proportional to the apparent degeneracy N . Finer details of the energy correlation inside the band are 
manifested for larger values of r. To understand the behavior of the spectral form factor, one writes the levels 
in the band (3 as ep(q) , q = 1, ■ • • , N, and substitutes in (0). Neglecting the cross-band correlations one gets 



N 



(7) 



This is the spectral form factor for a band, averaged over all the bands. The q summation can be performed by 
the saddle-point (or the uniform) approximation. The main contribution comes from the vicinity of the band 
extrema which correspond to the energy values where the spectral density is singular. That is, the prominent 
features in the form factor are due to the Van Hove singularities. Denoting by d^ep the second derivative of the 
band function at its extrema, one gets 

if(r) = C<(5 g %)- 1 ) /3 r- 1 , (8) 

where C is a numerical constant. It was shown in that the values of the constants which appear in (^J) and 
in (||) are compatible so that the two expressions match at r = 1/N. 

• t > 1 : The time interval is sufficient to resolve the point-like character of the spectrum. Hence, 

K(t)=j. (9) 



In the following sections we shall study how the expressions ( ||,||J9| ) make the transition to the Poisson form 
factor K(t) — 1 as disorder is introduced. The semiclassical (diagonal) approximation will be the starting point 
for the discussion of the transition in the first domain. This will be done in section II. To investigate K(t) in the 
second and the third domains, it suffices to study a system which has a single band in the periodic limit. The iV— site 
periodic Anderson model is such a system, and it will be discussed in section III. The important observation made in 
this section is that the transition is well described by considering the disorder perturbatively. The resulting explicit 
formulae for K(t) in the transition regime, reproduce the numerical data extremely well. The perturbative treatment 
also sheds light on the peculiar mechanism which reduces the value of K(t) from 7 to 1 in the third domain when the 
disorder splits up the degeneracies of the spectrum. We shall compare the results obtained separately for the three 
domains with numerical data for billiard and graph (network) systems. This will be done in section IV, where we 
shall summarize and discuss our findings. 



II. INTRODUCING DISORDER - THE SEMICLASSICAL APPROXIMATION 



We shall compute the spectral form factor (0) in terms of the Fourier transform of the oscillatory part of the spectral 
density. Using Gutzwiller's trace formula, d(r) can be expressed semiclassically as a sum over the periodic orbits j of 
the system 

d{T) = Y,^ T {r-r 3 )T J A J e 1 ^ (10) 

i 
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with primitive period Tj ~ r. Aj denotes the weight of the orbit corresponding to its stability and includes the Maslov 



phase. Sj is the action of the orbit in units of h. 
of repetitions of primitive orbits to the sum 



K{t) = 




Following the standard approximation, we neglect the contribution 
The form factor is now given by a double sum over periodic orbits 



i(sj-s,,) 



(11) 



It is well known 0], that for short time r this sum can be restricted to the diagonal terms j = f. However, when due 
to a symmetry, the orbit appears in gj different but symmetry-related versions, the contribution of all the symmetry- 
conjugated orbits must be added coherently. In such cases, (O) reduces within the diagonal approximation to 



K ( T ) ~ 9j S At(t - Tj) Tj \Aj 



(12) 



In the case of an extended, nearly periodic system, the diagonal approximation is valid up to r = 1/N , the Heisenberg 
time of the unit cell |^| . 

^From very general arguments it is clear, that in a system whose phase-space decomposes into several equivalent 
subspaces related by (unitary as well as anti-unitary) symmetry, the mean degeneracy g is just the number of such 
subspaces. Thus, if time-reversal invariance is the only symmetry obeyed, phase-space points with opposite momenta 
are equivalent and consequently phase-space is partitioned in g^ = 2 subspaces. In our problem, phase-space is 
invariant under a symmetry group containing N elements and therefore g = Ng^. Using the sum rule for periodic 
orbits MM, the form factor is finally written as 



K(t) ps g T NTP(r) 



(13) 



which we introduced in the previous section (|^). The normalisation of the staying probability P(t) is such that 
P(t) = Q,/w{t) at a time r, where the classical flow covers ergodically the part uj(t) of the total energy-shell volume 
f2. In particular, P(r — > oo) = 1 for an ergodic system and P(r) = N in a system which is composed of N unconnected 
ergodic cells. A more precise definition can be found in p]-lq]. For the present purpose, we need only the following 
property [^j: the return probability for a system composed of chaotic unit cells is independent of the presence or 
absence of long-range spatial order. Thus, within the diagonal approximation, the only effect of the introduction of 
disorder is the destruction of the coherence between the contributions of orbits which were related by symmetry in 
the original periodic system. This implies that in the diagonal approximation for K(r), (see (|l3|)) g = NgT, is to be 
replaced by g = g^. 

In order to describe the transition from g = Ng^ to g — <?t as the spatial symmetry is broken, we go slightly 
beyond the diagonal approximation ( [12] ) in that we retain in Eq. ( pT| ) the off-diagonal contributions from all those 
orbits which are degenerate in the symmetric system 



K(t) « 3t^<Sat(t - t t )t t \ A r \ 



N 

e i(A« r)j -As r ^,) 



(14) 



r runs now over all groups of symmetry-related orbits, while j, j' label the N orbits within each group. Possible 
degeneracies due to time-reversal are not affected by breaking the spatial symmetry and are thus contained in the 
prefactor <?t- The disorder which breaks the symmetry has been assumed weak enough such that (i) the orbits within 
the Heisenberg time of the unit cell t t < 1/N are structurally stable, i. e. no (short) periodic orbits appear or 
disappear due to the disorder, and (ii) the disorder does not alter by much the stability amplitudes and the periods 
within a group r so that in the prefactor A r _j ps A r , T rj - w r r . The variation of the actions are of the same order, but 
they cannot be neglected because they are measured in units of h and therefore, the resulting changes in phase, As rj - 
should be taken into account. 

Comparing ( |l2] ) and ( O ) we see that Eq. ( [i"3| ) represents the form factor also in the case of a weakly broken spatial 
symmetry, if g is replaced by an effective degeneracy 



9(r,S) 




(15) 
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which depends on the time r since the average on the r. h. s. is over all groups of periodic orbits r with length r r w r. 
The dimensionless parameter S has been introduced to characterize the strength of the symmetry-breaking disorder 
in a way to be specified in Eq. ( |l6| ) below. 

In order to evaluate Eq. (|15|), we need some information about the distribution of the disorder contributions to the 
phases of the periodic orbits As rj -. We assume, that the correlation length of the disorder is negligibly small compared 
to the mean length of an orbit. In this case As r j is a sum of many independent contributions, and the number of 
these contributions is proportional to the period of the orbit r. Hence, according to the central limit theorem, As r j 
are independent Gaussian random variables with mean value (As) = and variance 

(A 2 s)=S 2 t, (16) 

where the average is over all orbits of period r. With these assumptions we find from (plf ) 

g( T ,5) = ^( K N + N(N-l)\( C ^)\ 2 ) 

= g T (l + [N-l}c' 52T ). (17) 

In the first line we have used the fact that the N 3> 1 together with the statistical independence of As r j and As r ji 
justified above to replace the sum over j,f by its averaged value. In the second line, the Gaussian distribution of As 
was employed to give (e lAs ) = e _ ^ A s ^ 2 . It is easy to see, that Eq. ([l5]) indeed interpolates between g — Ng^ and 
g = <7t a s a function of the disorder. Note that the parameter which characterizes the disorder, S 2 , is multiplied by 
the time r over which the disorder acts. Hence, the classification of the disorder as "weak" or "strong" depends on 
the relevant time scale. 
In summary, we get, 

K(t, S) = g T (l + [N — 1] e- s2r )r P(r) = g T (l + [N - 1] e" 5 ^) V 'Nt /2 Cl ; t<1/N. (18) 

This expression provides the smooth transition from the periodic case, via the weakly disordered to the "metallic" 
domain. In section IV we shall show that this simple formula reproduces the form factor in the transition from 
periodicity to disorder very well. We emphasize once again that the present theory does not describe strongly 
disordered systems where the localization length is shorter than the system size. Such systems are outside of the 
scope of the present approach which is based on the "diagonal" approximation. 



III. INTRODUCING DISORDER - PERTURBATION OF A MODEL WITH A SINGLE BAND 

As explained in the introduction, the form factor in the domain r > 1/N is sensitive to the correlations among the 
levels which belong to the same band. Therefore, in order to investigate the form factor in this region, it is sufficient 
to study a model with a single band, which is what we do in the present section. In order to use the results of this 
section in the general context, we have to remember that the form factor in realistic systems is obtained as an average 
over many bands (see (Q)). This will smooth out several features of the single- band form factor, as will be explained 
in the sequel. 

The system we consider is a chain of iV unit cells of length a = 1, with periodic boundary conditions at the end 
of the chain. The chaotic scattering process in each cell is represented by a random potential and the dynamics is 
discretised on a lattice. Choosing convenient units, the Schrodinger equation reads 

- (<^n+l - 2^>„ + 4> n -l) + V n <j> n = E(j> n , (j) n = (j) L+n , (19) 

where 4>n is the wave function on the nth site. The on-site potentials V n are uncorrelated, random variables which 
are picked out from the same Gaussian distribution function with variance a. They obey 

(V n ) = 0, (V n V m ) = 5 nm a 2 n,m = Q...N-l, (20) 

The complexity of the scattering process is incorporated by neglecting the correlations between the potentials on 
different sites. 

In the periodic limit (cr = 0) , the levels are arranged in a discrete band 

E\ = 2(1 - cos(27rg/A0) 9 = 0, N - 1. (21) 
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The level density (compare Eq. [I]) 

w-i 

d(E) = J2 6 ( E ~ E J 

q=0 

exhibits van Hove singularities at the band edges E — and E = 4. This is a direct consequence of the periodicity of 
the system. 

For a ^ the singularity is smoothed out and for large values of a the level density becomes uniform between the 
upper and lower ends of the spectrum. This is the typical behavior we expect in generic one-dimensional disordered 
systems. ^From here on we consider the periodic case as the limit of the disordered system when a — ► 0. Accordingly, 
levels can be unfolded with the constant density. Since we consider here only weak disorder the mean level density is 
taken as 



(d) PS N/4. 



(22) 



The spectral form factor was defined in ([|). In the periodic case the energies are given by (|2l|). The spectral form 
factor is 



K(r,a = 0) = ± 



N-l 

exp(~iirE°rN/2) 

q=0 



(23) 



The second argument of the form factor denotes the strength of the disorder, and in the present periodic case it is 0. 
The expression (|23|) can be rewritten by expanding the exponential into a Bessel series: 



iNr-K cos -jr<} 



(24) 



fe=-c 



Exchanging the order of the summation over q and k yields 



K(t,0) = N 



OO 

J (nTN) + 2j2i nN JnN(nTN) 



(25) 



The function K (r, 0) which is shown in Fig. Q displays different features in the three domains of r. In the domain 
< t < 1/N the first term in ( p5| ) is dominant and Jq{tttN) sa 1. Hence, the form factor assumes the constant value 
N, and does not show any structure at all because we are dealing with a model with a single band. At r > 1/N, the 
form factor is a highly irregular function. It fluctuates more rapidly with increasing number of sites N. However, in 
order to compare the present theory with results which are derived for realistic systems, one should remember that in 
the latter case, the form factor is averaged over many bands which differ in their widths and structure. Such averaging 
can be effectively achieved by smoothing K(t, 0) over a small r window. 

The smoothed form factor (K(t, 0)) t is shown in Fig. [?]. In the range 1/N < r < 1/n ( p5|) is dominated by the 
Bessel function with zero index. The average behavior for large N and r < l/n can be approximated as 



(if(r,0)) T « iV(|J (^TAf ) r « 

7T Z T 



where we used the asymptotic form of the Bessel function 



(26) 



cos(z — vtt/2 — tt/4) 



and the average (cos 2 ) T = 1/2. 

In the third domain, r > l/n, the window-averaged function (K(r, 0)) T converges to a constant value. This constant 
is 7, the average degeneracy of the levels, and it approximately equals 2 since most levels are doubly degenerate (except 



Eq — and also E 



N/2 



4 if TV is even). In this range of r values, all Bessel functions contribute. Resumming the 



asymptotic forms of the Bessel functions we get 



(K(t,0)) t 



(27) 



G 



We conclude from this discussion that the one-band model, after proper averaging, reproduces the expected features 
of K(t, 0) in the relevant range r > 1/N. 

Introducing disorder, the form factor is given by 



K(r, a) 



N-l 



exp(-iirE°TN/2) 

9=0 



(28) 



where (. . ) a represents the average over the disorder and E° are the eigenenergies of the Eq. (y_9|) with cr -C 1. 

In the case of weak disorder, we can use degenerate perturbation theory to calculate how doubly degenerate energy 
levels are split. In first order the eigenenergies are given by 



e: 



±q 



Y V n exp(±i2—qn) 

n=0 



for q ^ 0, N/2 



E° for q = (and q=N/2 if N is even). 



(29) 
(30) 



Here we ignored a ^-independent constant, since it does not affect the form factor. The main effect of the perturbation 
is that it breaks the degeneracy of the energy levels which are symmetrically placed about the center of the band. 
The change of the mean level spacing is small, and can be neglected to leading order. 

Substituting the perturbed energy levels into ( p8| ) and leaving out the unimportant q = (and q = N/2 if N is 
even) levels we have 



K(r,a) 



1 

N 



N-l 

J2 2e~™ E l TN / 2 i 

q=l 



N-l 
n=0 



The disorder averaging can be performed analytically, as described in the appendix. For large N it leads to 

K(t, a) = l + A(2a) + A 2 {o){K{t, 0) - 2), 



(31) 



(32) 



where the universal function A{a) is defined in the appendix. A new combination of the variables involving the 
disorder strength shows up in this expression 



a 



7TT(7"\/ N 



(33) 



governing the properties of the transition from the periodic to the disordered case. For large r values the form factor 
converges to 1, since the perturbation breaks the degeneracy of the levels of the periodic system. Please note that- 
as in the semiclassical result Eq. ([i7|)- the deviation from the periodic form factor is governed by a dimensionless 
parameter containg the product of disorder strength and time. 
Approximating K(t, 0) by its average (p6h yields 



(K(r,a)) T 



(1 - e~ a y + ^ 
1 + A(2a) 



for r < 1/tt 
for t > 1/tt 



(34) 



The r < 1/tt part describes how the band structure is destroyed while the r > 1/tt part describes how the double 
degeneracy of levels is resolved (see Fig|8|). 

We can interpret the result for r > 1/tt in terms of the distribution p a (s) of splittings of levels s — E° — E°_ q which 
are degenerate in the periodic case. For large r the form factor is the Fourier transform of this distribution: 



K(j,a) = l+ f 
Jo 



dsp a (s) cos(stttN/2). 



Using the derived expression ( p6| ) for A(a), we can conclude that the splitting distribution has the form 

p a {s) = S -pw(sS), 



(35) 



(36) 



where pw{s) — nse 7rs / 4 /2 is the Wigner surmise and S = cry/~N is the mean splitting of levels. The Wigner surmise 
is known to be the exact distribution of the difference of the two eigenvalues of a 2 x 2 GOE random matrices. We 
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can conclude that in the present case, the ensemble of 2 x 2 matrix describing the splitting of levels in the first- 
order degenerate perturbation theory, reproduces the spacing distribution of the corresponding GOE with mean level 
spacing S. 

To check the applicability of the leading order perturbation theory, we computed the form factor numerically 
and compared with the analytical result. The parameter range was a = 0.002 ... 0.256 and N — 32... 256. The 



numerical results has been averaged for 1000 different disorder realizations. In Fig. |8| we compare formula ( |32[) and 
the simulations. We have found surprisingly good agreement in the whole range of r. The fact that K(t) displays a 
minimum where its value is less than 1, and that it approaches 1 asymptotically from below, is a direct consequence 
of the Wigner distribution of level splittings. In the next section we shall show that this formula applies very well 
also in the case of a multi-banded spectrum, indicating that the splitting distribution follows the Wigner distribution 
in more complicated situations too. 

IV. COMPARISON WITH NUMERICAL RESULTS AND DISCUSSION 

A. A chain of chaotic billiards 

The first class of systems which were investigated numerically are chains of chaotic billiards, (see Fig. |]) which can 
be arranged in a periodic (Fig. |l|(a)) or a disordered (Fig. |l|(b)) fashion. We denote the size of an individual billiard 
(i.e., the unit cell in the periodic case) by a and the chain length by L (and N = L/a 3> 1). In the following we discuss 

weakly disordered chains and assume that the conductance of the chain c = Nci ~ 1 . On time scales larger than the 
classical ergodic time for a single cell, the classical dynamics in the chain of billiards is diffusive, characterized by the 
diffusion constant D. In the diffusive regime, the classical dynamics of the system, and hence the diffusion constant, 
are, to a good approximation, independent of the strength of the disorder. The correlations in the quantum spectrum, 
however, crucially depend on whether the system is periodic or not, as discussed above. 

In this section we present numerical results for the form factor K(r) in weakly disordered chains. The case of 
periodic chains was analyzed in detail in two previous publications [pj. Here we focus on the crossover from the 
periodic case to the weakly disordered (metallic) case, which is predicted to follow (|l^) as a function of the disorder 
parameter 5. Due to time-reversal invariance of the billiard chain, <?t = 2. 

We have considered a chain composed of unit cells as shown schematically in Fig. ||. The sizes of the half disks were 
chosen so that the contribution of direct trajectories to the conductance is minimized. Disorder was introduced by 
shifting the disks at random to the right or to the left by a small amounts Ax. The dimensionless variance ((Ax) 2 ) 
is a measure of the disorder strength. 

Consider a periodic orbit j which hits M disks, m — 1, . . . , M. Its action Sj (measured in units of K) is affected by 
the M shifts Ax m and changes by an amount Asj, It is plausible that (A 2 Sj) cx ((kAx) 2 ) t and thus 

S 2 = C d ((kAx) 2 ) . (37) 

We have used (|37]) to estimate the quantity S 2 in ( |l8|) . The constant of proportionality Cd in ( |37| ) remains undeter- 
mined, it depends on the geometry of the system and on k. 

We have performed quantum-mechanical calculations for systems composed of N = 16 unit cells, with disorder 
parameters covering the domain of applicability of ([l8]) . The quantum-mechanical wave functions satisfy the Helmholtz 
equation augmented with Dirichlct boundary conditions on the channel walls and periodic boundary conditions along 
the chain. The quantum spectrum of this system can be determined using the method described in jl|] . In this way 
we have obtained the quantum spectra for several realizations of disorder, as well as for the periodic chain. 

Fig. ^ summarizes the results of our numerical calculations. It shows K(t, 5) as a function of r, in the periodic 
case, for weakly broken periodicity (four different disorder strengths) and for weak disorder. The conductivity per 
unit cell is independent of the disorder, and its numerical value was determined from a simulation of the classical 
dynamics of the system. A fit to the diffusion propagator at times larger than the ergodicity time allows one to 
determine D from which c\ ~ 33 emerges. The calculations were conducted for 6 values of the disorder strengths 
A 2 = (k 2 (Ax) 2 ) = 0.0, 6.87 x 10~ 4 , 9.55 x 10~ 4 , 1.49 x 10~ 3 , 2.56 x 10~ 3 and 4.38 x 10~ 3 . In all cases, we have 
calculated K (t) from 1500 eigenvalues, in the domain of k values which support 28 open transverse channels. 

The semiclassical theory for the periodic case reproduces the numerical results uniformly well over the three ranges 
of r values. The semiclassical theory matches very well the numerical results for the disordered systems in the domain 
r < 1/N. However, the numerical results in the domain l/N < r < 1 are not sufficiently smooth to allow a meaningful 



comparison with the theory developed in section (III). In this domain, the spectrum is afflicted by frequent near- 



degeneracies which make the calculation rather costly in terms of computer resources. This problem is circumvented 
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in the periodic case, where the translational invariance is used to facilitate the calculations. For larger values of the 
disorder, the degeneracy disappears, but, the effect we are interested in disappears, too. In the next subsection we 
discuss a different model exhibiting a transition from periodicity to weak disorder, where a quantitative comparison 
in the transition regime is possible. 



B. A chain of quantum graphs 

In this section we investigate a second model system — quantized graphs which were recently shown to provide 



an excellent example for a quantum chaotic system |10|. The graphs are defined by v = 1,...,V vertices and 



b = 1,...,B bonds with lengths L b connecting them. The wave function on a graph is a £?-component function 
(tpi(xi), • • • , iPb(xb)) ■ Each component satisfies the Schrodinger equation (ft = 2m = 1) 



d( , , r ) i> b {x b )=0. (38) 



At the vertices, the wave function must satisfy boundary conditions which impose continuity and current conservation. 
They guarantee that the Schrodinger operator is self adjoint, and its spectrum consists of discrete points. Implementing 
the boundary conditions, one derives a secular equation which provides a convenient means to compute the spectrum 
numerically. The graph is essentially a one-dimensional system, and therefore, the mean (wave number) spectral 
density is constant, proportional to the total length of the graph. 

The graph representing one unit cell was chosen to be the "cylindrical" network shown in the inset of Fig. |^. The 
cylinder consists here of n x — 2 layers with n y = 4 vertices each. The unit cell was constructed from more than one 
layer in order to remove any residual symmetry. Two bonds lead from each vertex to the neighboring layers, two 
more to other vertices in the same layer. Hence we have for the unit cell V = n x n y and B — 2n x n y . The lengths of 
all bonds are random, but the total length of the graph was fixed at Lr- = 2w such that the mean length of a bond 
is L = n(2N n x n y )~ l and the mean level spacing with respect to the wave number k is unity. For this reason it is 
natural to use, instead of energy and time, the wave number k and the length I as conjugate variables, since then no 
unfolding is necessary. In complete analogy to (|^) we introduce the spectral form factor via the length spectrum of 
the oscillating spectral density d(k) — Y] q 8{k — k q ) — 1 using a rectangular window 

K(r) = ^{\d(r)\ 2 ), (39) 



rk+Ak/2 

d(r) = / dfc'e- 27rifcr d(fc) . (40) 

t = l/2ir is simply given by the path length I measured in units of the Heisenberg length Lh = 27r. 

The classical analogue for the quantum graph is the random walk of a particle moving freely along the bonds and 



scattering at the vertices according to the quantum transition probabilities |10|. In the graphs we consider here, 
exactly four bonds are attached to each of the vertices. In this case the transition probability is 1/4 for all bonds, 
and the Lyapunov exponent is In 4 when time is scaled with the mean time between successive vertex traversals. The 
coarse-grained classical evolution is diffusive (n^) = D n n — Dil, where n w is the distance along the chain measured 
in unit cells (i. e. n w is the winding number in the periodic case), I is the length of a trajectory and n = l/L. 
When allowance is made for the fact that only half of the traversed bonds contribute to the diffusive transport, the 
diffusion constant is easily found from the analogy to a random walk on a ID-lattice with discrete time: D n = l/2n x , 
Di = N n y /irn x . 

The return probability entering (|l^) decays as r" 1 / 2 until it saturates at 1 when the diffusion covers the whole 
chain ergodically. The number of unit cells N = 8 was chosen such that this saturation occurs beyond the Heisenberg 
time of the unit cell t^ uc ' ) — 1/N and is thus not relevant for the form factor. In this case the return probability is 
explicitly given by 

P(r) = -4= (r c ( " c) < r < 1/N) (41) 



2-k^Dit 

Using g — 2N for the mean degeneracy of periodic orbits we finally obtain for the form factor 
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/ Ti 

K(t)=N,—Nt (Nt<1), (42) 

V TTTly 

which is shown in Fig. || with a smooth solid curve and has to be compared to the data obtained numerically without 
disorder (upper fluctuating curve). Beyond Nt = 1 the smooth curve shows the decay of the form factor as 1/r. 
Although the quantitative agreement is not perfect, the theory reproduces the essential features of the form factor, 
and in particular the peak at the Heisenberg time is correctly predicted. 
The disorder was introduced by small changes in the lengths of all bonds 

L <A) = L (o) + ALb (43) 

such that the total length remains constant (ALb) — 0. Here, b runs over all the NB bonds of the whole system. The 
strength of the disorder is characterized by the dimensionless parameter 

A 2 = k 2 (A 2 L b ) . (44) 

As shown in Fig. ^, the peak which characterizes K(t, 6 = 0) disappears gradually, when the strength 6 of the disorder 
is increased. In order to be able to apply the theory developed above, we have to take into account a feature which 
is particular to the graphs system. A periodic orbit of length r traverses on the average n(r) = tL^/L = ^Nrn x n y 
bonds, which, for sufficiently large r, can justify the discussion preceding ( ^6|) in section (|j). We have to bear in 
mind, however, that in fact not all of the n(r) length variations ALb accumulated in this way need to be independent, 
since in general, some of the bonds are traversed several times and moreover, for time-reversal symmetry the reversed 
bond contributes the same variation ALb — AL^. For this reason we introduce an average bond multiplicity m(r) for 
an orbit of period r. Then, the action variation of such an orbit is the sum of n(r)/m(r) independent contributions, 
each with a variance m 2 (r)A 2 . Hence we find for the variance of the sum 

(A 2 s) = n( T )m(T)A 2 . (45) 

In order to obtain an estimate for m(r) , we assume that a typical orbit covers ergodically some region of the phase- 
space such that each bond is traversed twice on the average (with momentum ±1) and hence m = 2. This is the case, 
e. g., at the Heisenberg time for an isolated unit cell, and — lacking a satisfying theory for m(r) — we have no choice 



but to generalize this special case. Comparing ( |45| ) with (16), we find for the disorder strength 6 2 = 8N n x n y A 2 . 
This is the parameter which we have chosen in Fig. ^ in order to compare the numerical data from Fig. |H| with the 
result of section O. In order to better distinguish the curves for small A we plot the quantity 1 — K§{t)/ Kq(t) which 
is according to ( O ) and ( [j"7|) given by (N — l)/N(l — e~ s ~ T ) and find indeed a reasonable agreement between the 
theory and the data. 

In Fig. [l^ we compare the graph data with the perturbative theory for a single band developed in section III. Since 
in our numerical calculations the number of unit cells N = 8 was not very large, we have to take into account the fact 
that for even N two levels in each band — at the border and in the center of the Brillouin zone — are not degenerate. 



Only the remaining N — 2 levels are described by the perturbative theory of section ( [ill]) , and consequently Eq. (|34| ) 
is replaced by 

K (t) = 1 + ^^A(2a) , (46) 

such that the asymptotic value in the periodic case is 2 — 2/N. Qualitatively, Eq. ( f46| ) predicts that the form factor 
for the periodic case has a minimum and beyond that approaches its asymptotic value from below. This non-trivial 
behavior is indeed observed in our numerical data. For a quantitative comparison we had to determine the unknown 
constant a which relates r to a according to Eq. (|33|). We have chosen a such that the position of the minimum 
in K(t) is the same for theory and numerics. Indeed this leads to a satisfactory agreement of Eq. (|i^) with the 
data, in particular beyond the minimum. It is reasonable, that this agreement becomes worse for smaller r, since 
then the main assumption behind Eq. ( [i"6|) — the lack of any correlation between different pairs of nearly degenerate 
levels — breaks gradually down. 

Summarizing our findings, we can confidently state that the numerical results displayed above provide convincing 
evidence in favor of the applicability of the simple semiclassical and perturbative approaches. This theory grasps the 
essential features of the transition, and provides simple expressions ([18|,[34]) for the form-factor and its dependence on 
the disorder. The main drawback of this theory is that it makes use of different approximations, depending on whether 
r is larger or smaller than the Heisenberg time 1/N. In the periodic limit, one could check the applicability of the 
theory in the vicinity of the Heisenberg time, by comparing it with the field theoretical expression which was derived 
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for periodic systems which violate time-reversal symmetry. The field-theoretical treatment |2| provides an expression 
which is uniformly valid for the entire r domain. The semiclassical theory of coincides with the field-theoretical 
expression in the separate domains of its validity, and did quite well even when the two expressions were extrapolated 
to the domain r w 1/N . A similar field-theoretical treatment of the the transition from the periodic to the disordered 
case does not exist yet, and it is naturally called for. 
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V. APPENDIX 



For calculating the averaged form factor in the disordered case, one needs the following quantities: 

(-l) k (x 2k ) 



(cos{x)) = ■ 



k=Q 



and 



(cos(x) cos(a;')) = 

fc=0,2=0 



(2*0! 



(-l) k+t (x 2k x' 21 } 
(2fc)!(2Z)! 



(47) 



(48) 



where 



N 



N ^ ( ^ \ 

14exp I i2— qn I 

71=0 ^ ' 



-N/2 



(49) 



and x' denotes the same, except the q is substituted by q' . Since the Taylor series of the cosine contains only 
even powers of its argument, after some simple but tedious calculations using the properties of Gaussian random 
distributions, we have the closed form: 



(X 2k ) = (TTT<j/2) 2k N k kl . 

The same calculations also show, that if q ^ q' then the variables are uncorrelated: 

(x 2k x' 21 ) = (7TTa/2) 2k+2l N k+l kW. = (x 2k )(x /21 ). 

Introduce the parameter 

tttv Act 



and using the property (J5l|), we define the function A(a) which appears in the expression ( |32| ) for K(t, ct): 

(cos(x)) = A(a). 

One can also show that 

(cos(a;) cos(n/)) = (cos(a:)) (cos(a: / )) = A 2 {a) 
(cos(x) cos(.t)) = / i + i C os(2x)\ = i(l + A(2a)) 



(50) 
(51) 

(52) 

(53) 

(54) 
(55) 



After substituting (|50|) into (E^) for the A(a) function results p5|: 
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M<») = E(- 1 ) fca2 '°7 2 lv = 1 - W° ? ' A Erfi(|o|/2) = 1 + V^-f e^ 4 Erf(ia/2) (56) 

where Erfi(x) denotes the error function for imaginary argument, the Erf(x) is the commonly used error function. 
The behavior of the function (see Fig. ||) for small arguments is Gaussian: 

A{a)^e a2 /\\-^ + o{a e )) (57) 
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(b) 



FIG. 1. Periodic (a) and aperiodic (b) chains of chaotic billiards. The chain length is denoted by L, a is the size of an 
individual billiard. Thus N = L/a is the number of units in the chain. 



FIG. 2. Typical discretized band spectra of a periodic chain with N = 16 unit cells. The energy levels are shown as a 
function of the Bloch phase 6 n for 10 bands in the case of (a) low, (b) intermediate and (c) high conductance. 
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FIG. 3. Unit cell of the chain of chaotic billiards. 
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FIG. 4. K (r) for chains of chaotic billiards and different strengths of disorder. 
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FIG. 5. The form factor for a quantum graph consisting of N=8 unit cells with 2x4 vertices each. The unit cell is shown 
in the inset with each dot corresponding to a vertex. The form factors were computed using the lowest 10,000 bands. Each 
spectral window in Eq. ( pj| ) contained 30 bands. The disorder strengths were (top to bottom) A = 0, 0.02, 0.05, 0.1, 0.5. The 
smooth curve represents the theoretical prediction without disorder. 
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FIG. 6. Reduction of the form factor due to disorder for various values of the disorder strength 

A = 0.01, 0.02, 0.05, 0.1, 0.2, 0.5 (bottom to top) compared to the prediction of Eq. @. 
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FIG. 7. The form factor (see Eq. 23) for the unperturbed system with lenght N = 64 (dots), the corresponding smoothed 
data (points) and the approximation Eq. 26 and 27 (lines). For r < 1/n the first term in the Bessel function expansion 
dominates the form factor. After averaged over 1/nN in r, the smoothed data fit to the theoretical function (|2(j). For r > 1/tt 
the soothed form factor converges to 7 ~ 2. Here the average is taken over a window of 1/n in r. 
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FIG. 8. Comparison between the form factor averaged for 1000 samples (a — 0.008, N = 32) and our perturbative expression. 
Both of them are smoothed over a 0.5 window in r. Points are from the numerical calculation, the solid line is the perturbativ 
formula of Eq. (Eh 
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FIG. 9. The universal function A(a) 
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FIG. 10. Large-r behaviour of the form factor for a periodic and a weakly disordered (A = 0.01) quantum graph (upper and 
lower solid lines, respectively) with N = 8 unit cells. The horizontal dashed lines represent the asymptotic values K (r) = 2—2/N 
and K (r) = 1. The heavy dashed lines show (K(r, er)) r according to Eq. (pi|). The parameter a entering Eq. ( |46[ ) via Eq. ( p3| ) 
has been determined by adjusting the the location of the minimum of the function A(2a) to the minimum observed in the 
numerical data. 
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